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Besides the conventional transverse couplings between superconducting qubits (SQs) and electromagnetic 
fields, there are additional longitudinal couplings when the inversion symmetry of the potential energies of the 
SQs is broken. We study nonclassical-state generation in a SQ which is driven by a classical field and coupled 
to a single-mode microwave field. We find that the classical field can induce transitions between two energy 
levels of the SQs, which either generate or annihilate, in a controllable way, different photon numbers of the 
cavity field. The effective Hamiltonians of these classical-field-assisted multiphoton processes of the single¬ 
mode cavity field are very similar to those for cold ions, confined to a coaxial RF-ion trap and driven by a 
classical field. We show that arbitrary superpositions of Fock states can be more efficiently generated using 
these controllable multiphoton transitions, in contrast to the single-photon resonant transition when there is only 
a SQ-field transverse coupling. The experimental feasibility for different SQs is also discussed. 
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I. INTRODUCTION 

Superconducting qubit (SQ) circuits M possess discrete 
energy levels and can behave as artificial “atoms”. In con¬ 
trast to natural atoms, with a well-defined inversion symme¬ 
try of the potential energy, these artificial atoms can be con¬ 
trolled by externally-applied parameters (e.g., voltage or mag¬ 
netic fluxes) and thus the potential energies for these 

qubits can be tuned or changed from a well-defined inversion 
symmetry to a broken one. Artificial atoms with broken sym¬ 
metry have some new features which do not exist in natural 
atoms. For example, phase qubits do not have an optimal 
point II0, so for these the inversion symmetry is always 
broken. 

When the inversion symmetry of these artificial atoms is 
broken, then the selection rules do not apply EMI, and 
microwave-induced transitions between any two energy lev¬ 
els in multi-level SQ circuits are possible. Thus, multi¬ 
photon and single-photon processes (or many different pho¬ 
ton processes) can coexist for such artificial multi-level sys¬ 
tems 03 , 111 111 - Two-level natural atoms have only a trans¬ 
verse coupling between these two levels and electromagnetic 
fields. However, it has been shown El that there are both 
transverse and longitudinal couplings between SQs and ap¬ 
plied magnetic fields when the inversion symmetry of the po¬ 
tential energy of the SQ is broken. Therefore, the Jaynes- 
Cumming model is not suitable to describe the SQ-field inter¬ 
action when the inversion symmetry is broken. 

Recently, studies of SQ circuits have achieved significant 
progress. The interaction between SQ circuits and the elec¬ 
tromagnetic field makes it possible to conduct experiments of 
quantum optics and atomic physics on a chip. For instance, 
dressed SQ states (e.g., in Refs. Eld) have been experi- 
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mentally demonstrated [18, 19]] . Electromagnetically-induced 
transparency (e.g.. Refs. |2Cj-(26|]) in superconducitng systems 
has also been theoretically studied. Moreover, Autler-Townes 
splitting EMI and coherent population trapping |[33l l have 
been experimentally demonstrated in different types of SQs 
with three energy levels. Experiments have shown that SQs 
can be cooled (e.g., Refs. 1134143711 ) using similar techniques 
as for cooling atoms. Moreover, sideband excitations H38LI41I1 
have been observed experimentally [421 l43ll using supercon¬ 
ducting circuits. Thus, SQs can be manipulated as trapped 
ions (e.g., in Ref. fH-Ell), but compared to trapped ions, the 
“vibration mode” for SQs is provided by an LC circuit or a 
cavity field. 

In trapped ions JHEl, multi-phonon transitions can be re¬ 
alized with a la_ser field. Multi-photon processes in SQs with 
driving fields [47] have been experimentally observed (e.g., in 
Refs. 1148145 311 1 when the inversion symmetry is broken. Thus, 
here we will show how nonclassical photon states can be gen¬ 
erated, via multi-photon transitions of a single-mode electro¬ 
magnetic field in a driven SQ, when the longitudinal coupling 
field is introduced. We will derive an effective Hamiltonian 
which is similar to the one for trapped ions. The single-mode 
quantized field can be provided by either a transmission line 
resonator (eq^., Refs. 091 |40L l54l l55ll T or an LC circuit (e.g.. 
Refs. |.V7. 5ii|). where the SQ and the single-mode field have 
both transverse and longitudinal couplings. In contrast to the 
generation of non-classical photon states using a SQ inside 
a microcavity lf58l46lll with only a single-photon transition, 
we will show that the Hamiltonian derived here can be used 
to more efficiently produce nonclassical photon states of the 
microwave cavity field when longitudinal-coupling-induced 
multiphoton transitions are employed. 

Our paper is organized as follows. In Sec. II, we derive an 
effective Hamiltonian which is similar to the one for trapped 
ions. We also describe the analogies and differences between 
these two types of Hamiltonians. In Sec. Ill, we show how 
to engineer nonclassical photon states using the multi-photon 
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FIG. 1: (Color online) Schematic diagram showing how two energy 
levels change with the external parameter for superconducting qubits. 
Here Eo and E\ are the eigenvalues of the ground and excited states, 
respectively. These vary with external parameters. For charge and 
flux qubits, the external parameters are the electric voltage and mag¬ 
netic flux, respectively. At the degenerate (or optimal) point, where 
the external parameter takes a particular value, the energy splitting 
reaches a minimum, A = hoj x , where the double potential well is 
symmetric. In this case, there is only a SQ-field transverse coupling. 
However, when the external parameter deviates from this point, the 
double potential well is asymmetric, and there are both transverse 
and longitudinal couplings between the SQs and the applied electro¬ 
magnetic field. 

coupling between the driven SQ and the quantized field. In 
Sec. VI, we discuss possible experimental implementations of 
these proposals for different types of SQs. Finally, we present 
some discussions and a summary. 

II. MULTI-PHOTON PROCESS INDUCED BY A 
LONGITUDINAL COUPLING 

A. Theoretical model 

As schematically shown in Fig.[Q the shape of the potential 
energy for some kinds of SQs (e.g., charge and flux qubits) can 
be adjusted (from symmetric to asymmetric and vice versa) 
by an external parameter, and thus the two energy levels of 
SQs can also be controlled. For charge and flux qubits, the 
external parameters are the voltage and the magnetic flux, re¬ 
spectively. However the potential energy of the phase qubits 
is always broken, no matter how the external field is changed. 
The generic Hamiltonians for different types of SQs can be 
written as 

h h 

H q = -oj z c r z + -uj x a x . ( 1 ) 

As in experiments, we assume that both parameters co z and 
uj x can be controlled by external parameters. The parame¬ 
ter uj z = 0 corresponds to the optimal point and well-defined 


inversion symmetry of the potential energy of the SQs. How¬ 
ever, both nonzero parameters and c o x correspond to a bro¬ 
ken inversion symmetry of the potential energy of the SQs. 
Below, we first provide a general discussion based on the qubit 
Hamiltonian in Eq. CD, and then we will specify our discus¬ 
sions to different types of SQs. The discussion of their exper¬ 
imental feasibilities will be presented after the general theory. 

Let us now assume that a SQ is coupled to a single-mode 
cavity field and is driven by a classical field, where the Hamil¬ 
tonian of the driven superconducting qubits is 

H = H q +huja^a+hga z (a+a^)+hfldcr z cos(LOdt+<f>d)- (2) 

Here, a'(a) is the creation (annihilation) operator of a single¬ 
mode cavity field with frequency w. The parameter O,/ is the 
coupling constant between the SQ and the classical driving 
field with frequency a;,/. The parameter g is the coupling con¬ 
stant between the SQ and the single-mode cavity field. The 
parameter <fi d is the initial phase of the classical driving field. 

Equation 0 shows that there are transverse and longitudi¬ 
nal couplings between the SQs and the electromagnetic field. 
This can become clearer if we rewrite the Hamiltonian in 
Eq. (0 in the qubit basis, that is, 

H = —uj q a z + huja^a + Hg z a z (a + a^) + hg x a x {a + a^) 
+ fr£l dz o z cos (u d t + </> d ) + h^dxcrx cos(w d f + (3) 

with four parameters g z = g cos 9, g x = —g sin 9, Q dz = 
Q d cos9, and n dx = —fid sin 9. Here, the parameter 9 is 
given by 9 = arctan (uj x /uj x ), and the qubit eigenfrequency 
is u q = ^/ujf+ui f. 

The Hamiltonian in Eq. ([3]) shows that the qubit has both 
transverse and longitudinal couplings to the cavity (driving) 
fields with transverse g x (fl x ) and longitudinal g z (fL) cou¬ 
pling strengths. When both uj z = 0 and Q d = 0, Eq. (0 is 
reduced to 

H = | u} q a z + Hga x (a + a f ) (4) 

which has only the transverse coupling between the SQ and 
the single-mode field. If we further make the rotating-wave 
approximation, then Eq. 0 can be reduced to the Jaynes- 
Cumming model, which has been extensively studied in quan¬ 
tum optics That is, there is only a single-photon transi¬ 
tion process when the qubit is at the optimal point. However, 
the transverse and the longitudinal couplings between the SQ 
and the single-mode field coexist, when the inversion sym¬ 
metry of the potential energy is broken and oj z is nonzero for 
the SQs. As shown below, this coexistence can induce multi¬ 
photon transitions between energy levels of SQs and make it 
easy to prepare arbitrary nonclassical states of the cavity field. 

Below, we assume that both u) z and u> x are nonzero. We 
also assume that the SQ and the quantized field satisfy the 
large-detuning condition, that is. 


u: q = i/u > x +wl » w. (5) 

In this case, the SQ and the quantized field are nearly de¬ 
coupled from each other when the classical driving field is 
applied to the SQs. 
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B. Multi-photon processes and sideband excitations with 


Let us now study how multi-photon processes can be in¬ 
duced via a longitudinal coupling by first applying a displace¬ 
ment operator 


D 




to Eq. © with 


(6) 


H 0 = -u z a z + huja^a, (14) 

we can further expand the Hamiltonian in Eq. CD, in the in¬ 
teraction picture, into 

H int = ^oj x J2 {JN n (t)<7+a^ m a n + h.c.}, (15) 

N,m,n 


V = 2 g/w. 


( 7 ) with 


Thus // is the normalized qubit-cavity coupling. It is also 
known as the Lamb-Dicke parameter. Hereafter, we denote 
the picture after the transformation D (j]a z / 2) as the displace¬ 
ment picture. In this case, we have an effective Hamiltonian 

H e s = DHD t = —w z a z + hujc^a + M}dcr z cos(u>dt + <f>d) 

+ {er+ exp [rj (a f - a)] + h.c.} . (8) 

From Eq. © with VL d = 0, we find that if nuj = w 2 , then the 
multiphoton processes, induced by the longitudinal coupling, 
can occur between two energy levels formed by the eigen¬ 
states of the operators a z . However, such process is not well 
controlled. Moreover, ui z is usually not perfectly equal to nco, 
for arbitrarily chosen n. These problems can be solved by 
applying a classical driving field, in this case Q d ^ 0. 

To understand how the classical field can assist the cavity 
field to realize multi-photon processes in a controllable way, 
let us now apply another time-dependent unitary transforma¬ 
tion 


U d (t) = exp 




(9) 


( 10 ) 


to Eq. © with the Hamiltonian H d defined as 

H d (t) = ^^ d<Tz sin (ui d t + 4> d ), 
u d 

and then we can obtain another effective Hamiltonian 

= U d H eS U ] d - iU= ^u z a z + fuva^a 

, OO 

H—— {JN<r+BN(t)+ H.c.}, (11) 

N=—oo 


where the time-dependent expression Bn (t) is given as 

B n (t) = exp [rj (al - a) + iN(u d t + 4> d )\ ; (12) 

./v = Jjv ('J'd) is the A’th Bessel function of the first kind, 
with x d = 2n d /uj d and J N {x d ) = (-1 ) N J- N (x d ), and ijis 
similar to the Lamb-Dicke parameter in trapped ions 1441451. 
Via the unitary 


V 0 (t) = exp 



(13) 


JTit) 


(-1 YJn +n 
mini 


exp 



(16) 


x exp [iN(u d t + 4> d ) + i (m — n) uit + iu> z t]. 


Equation (fl5l > clearly shows that the couplings between the 
SQs and the quantized cavity fields can be controlled via a 
classical field when they are in the large-detuning regime. 
Comparing the Hamiltonian in Eq. (IT5l > with that for the 
trapped ions 144 1451. we find that the Hamiltonian in Eq. (fl5l i 
is very similar to that of the two-level ion, confined in a 
coaxial-resonator-driven rf trap which provides a harmonic 
potential along the axes of the trap. Therefore, in analogy 
to the case of trapped ions, there are two controllable multi¬ 
photon processes (called red and blue sideband excitations, 
respectively) and one carrier process'. 


(i) when n > m, with n — m = k, and the transition satisfies 

the resonant condition Nui d = uj z — kto, with N, k = 
1, 2, 3, • • •, the driving frequency Noj d is red-detuned 
from the qubit frequency w z . Thus, we call this multi¬ 
photon process the red process. 

(ii) when n < m, with m — n = k, and the transition satisfies 

the resonant condition Nui d = uj z + kto, with N, k = 
1, 2, 3, • • •, the driving frequency Nuj d is blue-detuned 
from the qubit frequency uj z . Then we call this process 
the blue process. 

(iii) when n = m and lo z = Nuj d (N = 1, 2, ■ ■ ■), the 
driving field with N photons can resonantly excite the 
qubit. We call this transition the carrier process. 


However, there are also differences between the Hamilto¬ 
nian for trapped ions |44, j4*§] and that in Eq. ( | 1 5[ i. These 
differences are: 


(i) For a given frequency u> d of the driving field, there is 
only one multi-photon-transition process in the sys¬ 
tem of trapped ions to satisfy the resonant condition, 
but the SQs can possess several different multipho¬ 
ton processes, resulting from the longitudinal coupling 
between the classical field and the SQ. For instance, 
with the given frequencies u> d and u>, and for the cou¬ 
plings with the Nth and A''th Bessel functions, two 
transitions with the red sideband resonant conditions: 
Noj d = oj z — kuj and N'uj d = oj z — k'ui, might be sat¬ 
isfied. Once the condition (N — N')/(k' — k) = Lo/u> d 
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is satisfied, then these two resonant transitions can si¬ 
multaneously occur. Similarly, for the case of blue- 
sideband excitations, the condition that two resonant 
transitions simultaneously occur for the k and pho¬ 
ton processes is (N — N')/(k — k') = u)/cod■ We 
can represent the transition type in the sign of k and 
k'. Thus if we want some terms with N' unresonant, all 
we need to do is to let (N — N')/(k — k') ^ ui/uid, 
i.e., u>z 7^ ui ( N'k — Nk') / (N' — N) . One sufficient 
condition is that ui z ^ c on/ (N' — N) ( n = 0, ±1, ±2, 
•••)■ 

(ii) The Lamb-Dicke parameter rj for the trapped ions is 

determined by the frequency of the vibration phonon, 
mass of the ion, and the wave vector of the driving field. 
However the Lamb-Dicke parameter r/ here is deter¬ 
mined by the frequency ui of the single-mode quantized 
field and the coupling constant g between the single¬ 
mode field and the SQ. 

(iii) For multi-photon processes, the coupling between 
trapped ions and the phonon is always on. However, 
such processes can in principle be switched off at the 
zeros of the Bessel functions of the first kind. 

(iv) The term in Eq. dT5l > with N = 0 means that the driving 

field has no help for the excitation of the SQ. Thus this 
term is neglected in the following discussions. How¬ 
ever, the driving field can always be used to excite the 
trapped ions when certain resonant condition is satis¬ 
fied. 

(v) For trapped ions, the ratio between the transition fre¬ 

quency of the qubit and the frequency of the vibration 
quanta is often about 10 9 . Thus the upper bound for 
the photon number k in the multiphotn process is about 
k = 10 9 . However, in the SQ circuit, the frequency of 
the SQ can be several tens of GHz, and the quantized 
cavity field can be in the regime of GHz. Thus the pho¬ 
ton number k is not extremely large. For example, if 
uj z = 20 GHz and w = 2 GHz, then the upper bound 
for k is 10. 

To compare similarities and differences. Table [T] lists the 
main parameters of the Hamiltonian for trapped ions and those 
of the SQ in Eq. ( fTTb . We should note that the Lamb-Dicke 
parameter q can become very large in circuit QED systems 
in the ultrastrong | [7ll473ll and deep-strong lp74H77l l coupling 
regime. Our discussion below is in the ultrastrong coupling, 
but can be straightforwardly extended to the deep-strong cou¬ 
pling regime. 

C. Bessel functions and coupling strengths 

In the process of generating nonclassical photon states, the 
coupling strength J^ n {t) plays an important role. In our 
study here, the Bessel functions of the first kind are crucial 
factors in the coupling strengths. The possible values of the 


TABLE I: Comparison of some parameters between the Hamiltonian 
in Eq. 03 and that of the trapped ions (e.g., in Ref. [45j]). Here LD 
refers to the Lamb-Dicke parameter. 


Parameters 

Superconducting qubits 

Trapped ions (orders of 


(orders of magnitude) 

magnitude) 

LD parameters 

2 g/ui ~ (0.2 - 1.8) 

v ~ (0.2 - 0.9) 

Carrier Rabi frequencies 

Renormalized Jn u x / 2 

Renormalized Q 

Driving field frequencies 

u d (N = !,•••) 

UL 



FIG. 2: (Color online) Bessel functions Jn{xci) of the first kind, 
with N — 0, 1, 2, 3, are plotted as functions of the ratio Xd in (a). 

J ? 0 (t) |, | J 9 1 (f) |, and | Jj 1 (t) | have been plotted in (b), (c), and (d) 
as functions of g and Xd, respectively. Recall that Xd = 2 fld/oJd 
is the ratio between the driving field-SQ Rabi frequency fid and the 
frequency u)d of the driving field. 


Bessel functions depend on the ratio Xd between the driv¬ 
ing field-SQ Rabi frequency f Id and the frequency lo,i of the 
driving field. For several recent experiments with supercon¬ 
ducting quantum circuits, the coupling constant f Id is usu¬ 
ally in the range from several tens of MHz to several hun¬ 
dreds of MHz, e.g., 10 MHz < Q d < 500 MHz. The fre¬ 
quency LOd of the driving field is in the range of GHz, e.g., 
1 GHz < u: d < 20 GHz. Thus the ratio x d is in the range 

10~ 9 < x d < 1. (17) 

For completeness and to allow a comparison between them, 
several Bessel functions are plotted as a function of the pa¬ 
rameter Xd = 2f ld/uid in Fig. [2a), which clearly shows 
J 0 (x d ) > Ji{x d ) > h{x d ) > > JN^Xd ) in the range 

of 10 -9 < Xd < 1. Thus if the classical driving field is cho¬ 
sen such that the ratio x d is less than 0.5, then we only need 
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to consider the terms in the Hamiltonian in Eq. (IT5l) with the 
Bessel functions Jo(xd) and \ J 1 (x,i)\ = \ J-i(xd)\, and other 
terms are negligibly small. As discussed above, it should be 
noted that the frequency u>d of the driving field has no effect 
on the coupling between the SQ and the quantized field in 
terms of the Bessel function Jo(xd )■ Thus the driving-field- 
assisted transitions between the SQ and the quantized field are 
determined by the terms with the Bessel functions J±i(xd), 
when other high-order Bessel functions are neglected. Fig¬ 
ure [2ja) also shows that the terms with the Bessel function 
J± 2 {xd) are also not negligible when x,/ becomes larger, e.g., 
0.5 < x,i < 1. Thus in the regime 0 < x,i < 1. the terms with 
high-order Bessel functions (e.g., the ones with N > 3) can 
be neglected. 

As an example. Figs. [2b, c, d) illustrate how | are 

affected by to, n, Xd, and tj. Since the maximal point occurs 
at ij = \J tti T n, if other variables are fixed, thus we can find 
an obvious shift of the maximal point along the //-axis with 
increasing m + n. We can also find that \J^ n (t)\ have similar 
results as those of | J{ ran (i)| versus to, n, Xd, and ?/. By tun¬ 
ing Xd and rj, we can change the Rabi frequencies, and thus 
optimize the generation time. 


III. GENERATING NON-CLASSICAL PHOTON STATES 
USING SUPERCONDUCTING QUANTUM CIRCUITS 


In this section, we discuss how to generate non-classical 
photon states via transverse and longitudinal couplings be¬ 
tween SQs and the single-mode cavity field, with the assis¬ 
tance of a classical driving field. 


A. Interaction Hamiltonian and time-evolution operators 


Let us now analyze the interaction Hamiltonian and time 
evolutions for the three different processes based on the 
Hamiltonian in Eq. (1 1 5b . We have three different interaction 
Hamiltonians. In the interaction picture with the resonant con¬ 
ditions of different photon processes, by assuming n = m + k 
for the red-sideband excitation, m = n + k for the blue- 
sideband excitation, and m = n for the carrier process. We 
now discuss the general case for coupling constants with any 
number of Bessel functions. For a red process with the Nth 
Bessel functions, we derive the Hamiltonian 


Hr = £1 V ( ~ 1} 2, a+a .. a + h.c, (18) 


’N.r / j 


to!(to + fc)! 


with the resonant condition 


Nojd = w 2 — kco. 


For a blue process with the TVth Bessel functions, we have the 
Hamiltonian 


„ _ r (fc) V- (-l)V"o+a t( "+ fc )o" , , „ 

Hb ~ J N,b22 n \(n + k)\ +h ’ ( 9) 


with the resonant condition 

Nu>d = u> z + kui. 

The parameters J and J b r b for the red process in Eq. (ITHl 
and the blue one in Eq. ( fl9l > are given by 


r(fc) { r 

J iV,r = ( _1 ) ^- J ^ ex P 




T (k ) &x T 

J N,b = ex P 


1 f2g 


- iN<j^P 


( 20 ) 

( 21 ) 


where the subscript /3 takes either r or b, we use <j>^ to char¬ 
acterize the initial phase of either the red or the blue process. 
For the carrier process with the Nth Bessel functions, the in¬ 
teraction Hamiltonian is given by 


TT _ r(o) 

n c — j Nc 2_^ v 


(—1 ) n a + a^ n a n 


i\n\ 


+ h.c., 


( 22 ) 


with the resonant condition 

Nuj d = ijJ z , 

and the coupling constant 


t(0) 1 7 

j n,c = ~UxJn exp 


1 


--77 + iN(j> d 


(c) 


(23) 


We also note that all non-resonant terms have been neglected 
when Eqs. (118H22I) are derived. The dynamical evolutions of 
the systems corresponding to these three different processes 
can be described via time-evolution operators. For example, 
for the kth red, blue, and carrier sideband excitations, we re¬ 
spectively have the evolution operators 

k— 1 00 

U nI ( f ) = £ l n ) a °° + £ C0s ( ^ ( n \ an 


n —0 


n —0 


+ ^2 e l ^ N ' r l7r ^ 2 sin ^ tj | n + k) (n\ cr_ 

n=0 

00 

+ ^2 e^ J '’ r_l ’ r/2 sin ^ tj | n) (n + k\ a + 

n =0 
00 

+ ^2 C0S ( ^N n r *) \ n + k )( n + k \ &00, (24) 


n =0 
k -1 


U n\ (*) = £ ^ ( n l an + £ cos ( *) \ n ) ( n\ o-oo 


n =0 


n —0 


-£< 

n —0 
00 




sin ^ VL k ^ n b t^j | n + k) (n\ cr + 


+ £ e l< ^’ b l 7 r ^ 2 si n ( tj | n) (n + k\ tr_ 

n=0 
00 

+ ^2 cos ^ tj \n + k) (n + k\ an, (25) 


n —0 
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and 


i>( 

n°’ n 

“iV,c 

n =0 


oo 

+ E e_ 

Cc~ 


n—0 


+ 


+ 


oo 

£ eO^sin (|n^ c | t) | n) (n\ a+ 

n—0 

oo 

^2 cos ( tj \n) (n\ a 00 , (26) 


n—0 


where the complex Rabi frequency and its phase angle are 
respectively defined as 


^Pp = J N^P\j ^ ’ ^ 27 ' > 

<t>N*P = ar § (pPp) > ( 28 ) 

with /3 = r, 6, c and ,7^ are given in Eqs. (l20l >. QTk and 

(f23l> . Here iJ,!' 1 (x) represents the generalized Laguere poly¬ 
nomials. Let us assume that the two eigentates \g) and \e) of 
the Pauli operator a z satisfy cr z \g) = — \g), and a z \e) = |e), 
then we define the following operators an as (Too = | g) (<?|, 
o-u = |e) (e|, ctqi = \g) (e|, and cri 0 = \e) (g |. 


B. Synthesizing nonclassical photon states 

We find that the interaction Hamiltonians in Eqs. (1181122b in 
the displacement picture are very similar to those for trapped 
ions B45I1 . Therefore, in principle the non-classical photon 
states can be generated by alternatively using the above three 
different controllable processes. We expect that the prepared 
target state is 


22-max 

I</Vmax) =j2 Cn I”) 0 I 9 ) ’ (29) 

n—0 

where n max is a maximal photon number in the photon state 
of the target state. Here, | n) ® |(?) = | n) \q) denotes that the 
cavity field is in the photon number state | n) and the qubit is 
in the state | q), which can be either the ground | g) or excited 
|e) state. The parameter \C n \ 2 is the probability of the state 
| n) ® |g). The steps for producing the target state for both the 
case q = g and q = e are very similar. We thus take q = g as 
an example to present the detailed steps. The target state then 
takes the form 


22-max 

IV'nmaJ =J2 C k | k) ® I g) . (30) 

k -0 

We point out that all the states here (e.g., the target state) are 
observed in the displacement picture, if we do not specify this. 


We assume that the system is initially in the state |0) | g). 
Then, by taking similar steps as in Ref. 145jl, we can generate 
an arbitrary state in which the states \ijj n ) in the nth step and 
in the (n — l)th step have the following relation, 

l^n) = (t n ) ufy n (t„) (0) |^ n _i), (31) 


with 


\^n) = Y, C kg\k)\g)+C^me). (32) 
k =0 


Here if n = n max , C 0e = 0, and cj£ = C k , then \ip n ) in 
Eq. (l32l) is reduced to |V’n max ) in Eq. (f30t . Above, t n is the 
time duration of the control pulse for the nth step. The unitary 
transform (t) is defined as 

U^ n) it) = Ho it) U^ n) it) . (33) 

Here Vo it) is given in Eq. ( fl3l ). Also, UP it) is actually 
Ud it) in Eq. (|9]). but with uid and <pd replaced by coP an d 
, which denote respectively the frequency and phase of 
the driving field for the nth step. Moreover, (t n ) de¬ 

notes a unitary transform of the nth step, and is taken from 
one of Eq. (l24l >. Eq. (l25l) . and Eq. (l26l) depending on which 
one is chosen as /3 n among the characters “r”, “6”, and “c”. In 
(f„), the parameters uid and 4>d must also be replaced 

by and <2d \ respectively. 

The target of the nth step is to generate the state \ijj n ) from 
the state \ip n _i). We assume |-0 n _i) is in the displacement 
picture, which is the state generated after the (n — l)th step. 
We first use UP (t n ) to transfer IV’n-i) from the displaced 
picture into the interaction picture. Then we choose one of 
the evolution operators in Eqs. (f24l) (f26l) with a proper pho¬ 
ton number to reach the target state in Eq. ( BTI i. Since the state 
\ip n ) should also be represented in the displaced picture, after 
the state of the nth step via the evolution operators uipp (£ n ) 
and Up (0), we have to transfer it back to the displaced pic¬ 
ture, which results in the appearance of Up { (f n ) in Eq. ( |3T1 ). 

The longitudinal coupling results in multi-photon pro¬ 
cesses. Thus the state preparation using the longitudinal cou¬ 
pling is in principle more convenient than that using a single¬ 
photon transition in the usual Jaynes-Cumming model 115811 . 
For example, the Fock state | n) can be generated with a carrier 
process and a longitudinal coupling field-induced n-photon 
process. However, it needs 2 n steps (n step carrier and n 
step red-sideband processes) to produce a Fock state |n) if 
we use the Jaynes-Cumming model 1581. The selection of 
Upp (t n ) in Eq. (OH for each step is almost the same as 
that in Ref. ll45ll . That is, the target state in Eq. Ool ) can be 
obtained either by virtue of one carrier process and n max red- 
sideband excitations, or by virtue of one carrier process and 
n mi x blue-sideband excitations. 

The steps to generate the target state in Eq. O0l > from the 
initial state |0)|<?) using carrier and red sideband excitations 
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FIG. 3: (Color online) Schematic diagram for state generation of the 
target state |V>„ max ) = Z)fc=“ Ck l fc ) ® Ifl). for example, n max = 2. 
(a) Step 0, the system is initially in the state |0) ® |t/), and the arrow 
means that a carrier process is applied with the 0-photon inside the 
cavity, denoted by the operator (to), (b) Step 1, after the step 
0, the system is in the state Cq°J |0) ® \g) + Cq°J |0) ® |e), in which 
the parameters C^} and C ^ are determined by the time duration to- 
The arrow means that a 1-photon red process is applied with the time 
evolution operator ujj (ti) after the step 0. (c) Step 2, after the step 

1, the system is in the state |0)®|<?)+Ci^ |l)®|ff)+Cg^ |0)® 
|e), in which the coefficients of the superposition are determined by 
the time duration to and ti. The arrow means that a 2-photon red 
process Uffi (ti) is applied to the system after the step 1 . (d) After 
the step 2 (with well-chosen time durations to, t\ and ti), the system 
is in the state Co |0) ® |<?) + Ci |1) ® |g) + Ci |2) ® |g), which is 
just the target state \ip-z). Other superpositions can also be generated 
using similar steps. 


are schematically shown in Fig. fusing a simple example. All 
steps for the required unitary transformations are described as 
follows. First, the initial state |0) |g) is partially excited to 
| 0 ) |e) by a carrier process ( n = 0 ) with a time duration £q 
such that the probability ICg^j 2 in | 0 ) | < 7 ) satisfies the condi¬ 
tion |Cq ^| 2 = |Co| 2 , with Co given in Eq. (l30t . After the 
carrier process, the driving fields are sequentially applied to 
the qubit with n max different frequency matching conditions, 
such that a single-photon, two-photon, until 7i max -photon red 
processes can occur. Thus the subscript 3 rl in the unitary 
transform U^p satisfies the conditions j3 n = c with n = 0, 
and (3 n = r for n > 1. By choosing appropriate time dura¬ 
tions and the phases of the driving fields in each step, which 
in principle can be obtained using Eq. ( |3H . we can obtain the 
target state shown in Eq. ( l30l ). The detailed descriptions can 
be found in Appendix [A] 


IV. THE INITIAL STATE AND TARGET STATE 


Above, we assumed that the target state is generated from 
the initial state which is the vacuum state in the displacement 
picture defined by Eq. ([ 6 }. However, in experiments, the ini¬ 
tial state is usually the ground state, obtained by cooling the 
sample inside a dilution refrigerator. We now investigate the 
ground state of the effective Hamiltonian when there is no 
driving field. The Hamiltonian H e s without driving field can 
be expressed as 

H' eS = ^uj z a z +huja^a+^co x {0 + exp [77 (a* - a)] + h.c.} , 

(34) 

in the displacement picture. However, in the original picture, 
the corresponding Hamiltonian is 

H' = : u; z a z + ^ uj x g x + hwa)a + hga z (a + a t ), (35) 


which possesses the characteristics of broken-symmetry and 
strong coupling and is hence difficult to solve analytically. 
Due to the mathematical equivalence between Eq. (l34t and 
Eq. (f35l> . it is also difficult to solve Eq. d34l i analytically. We 
thus resort to numerical calculations to obtain the ground state 
of H' eS . We define the ground state of the Hamiltonian H' eS 
as \ip g ), and the probability of the ground state | i/j g ) to be in 
the vacuum state |0) as P g ,o- The relation between \ip g ) and 
Pgfi can be written as 


I i’g) = [& 10) + V 1 - IC 0 I 2 1 9 ), 
Pg ,0 = 1^01 2 , 


(36) 

(37) 


where 1 8i/) g ) denotes a superposition of photon number states 
except the vacuum state. In Fig. [4] as an example, by taking 
u> z /2tt = 19.5 GHz and 10 = 2 GHz, we have plotted P 9: 0 as 
a function of rj and oj x . We find that the probability P g o > 
0.99, at least in the region 0 < 77 < 3.5 and 0 < t o x /oj z < 
0.2. More specifically, the ground state of the Hamiltonian in 
Eq. ( l34l ) is closer to the vacuum state when the parameters 77 
and lj x are smaller. Thus, our assumption that the initial state 
of the cavity field in the displacement picture is the vacuum 
state, can always be valid only if the related parameters, such 
as oo x and 77 , are properly chosen. 

We have demonstrated how to generate an arbitrary super¬ 
position of different Fock sates from the vacuum state in the 
displacement picture. Thus, once the state is generated, we 
have to displace the generated state back to the original pic¬ 
ture via the displacement operator I) ' (r]a z /2). For example, 
the initial state | 0 ) | < 7 ) in the displacement picture becomes 


£> t (Jot) |0)|g) 



(38) 


in the original picture, where \a,n) = D (a) \n) denotes the 
displaced number state 0 . Similarly, the target state |t/V i max ) 
in the displacement picture becomes 


''■max 

KnJ = H |,n) Iff), (39) 


n —0 
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FIG. 4: (Color online) Probability P g ,o for the ground state to be 
the vacuum state in the displacement picture as a function of 77 and 
Here we assume uj z /2-k = 19.5 GHz, and ui = 2 GHz. Re¬ 
call that rj is the normalized coupling or the Lamb-Dicke parameter. 


in the original picture. It is obvious that the initial state 
\rj/2, 0) of the cavity field in the original picture is a coherent 
state with the average photon number ( 77 / 2 ) , while the target 
state is the superposition of the displaced number states. 

The statistical properties of a displaced number state with 
D(a)\n) can be described by the probabilities of the photon 
number distribution as below 


(l\a,n) = 


a 1 "v^TF r{l-n) (I ,2\ 

exp(MV2) L ” U a U’ 

l > Tl 

r(n-l) f 1 i2\ 

exp(|a| 2 /2) 1 V J 

, l < n 


(40) 


Thus the displaced target state in Eq. (I39l > can be written as 


,l max CXJ 

l<CJ= Z ^Z^x* II’ n )i ff > 

n =0 1=0 

00 

= E c * D i^>> ( 41 > 


i=0 


where C® = {l\r]/2, n). The probability of the 

target sate |'0^ , max ) to be in the photon number 1 1) in the orig¬ 
inal picture can be given as 


p, D = K \ 2 


n 


max ((max 

n =0 m =0 




()■ (42) 


In Fig. [5] as an example, the distribution probabilities Pf are 
plotted for different photon states, that is, |^^ max ) is taken 
as 177 / 2 , 0) | g) , 177 / 2 , 2 ) | g) , or (( 77 / 2 ,0) + ( 77 / 2 , 2 )) \g) /y/2, 
which is 10), |2) and (|0) + 12))/\Z2, respectively, in the dis¬ 
placement picture. Figure [5] shows that the photon number 
states in the displacement picture are redistributed after these 
states are sent back to the original picture. Even though the 


0.5 

0 

0.5 

or 

0 

0.5 


0 


FIG. 5: Photon number distributions of (a) 1 77 / 2 , 0), (b) 1 77 / 2 , 2), and 
(c) (|? 7 /2, 0) + 1 77 / 2 , 2)) /\/2. Here l refers to the photon number 
and Pi refers to the probability on |Z). We have taken the Lamb- 
Dicke parameter 77 = 2 g/Q. = 0.7. 



state in Fig.|5jc) is the linear sum of the states in Fig.[5ja) and 
Fig- 0b), the photon number distributions are not linearly ad¬ 
ditive. Because the interference between different displaced 
number states, which corresponds to the terms of m ^ n in 
Eq. (l42l >. can also give rise to the variation of the photon num¬ 
ber distribution. It is clear that a number state in the displaced 
picture can become a superposition of number states in the 
original picture, which might offer a convenient way to pre¬ 
pare nonclassical photon states. 


V. NUMERICAL ANALYSIS 

We have presented a detailed analysis on how to prepare 
nonclassical photon states using the longitudinal-coupling- 
induced multi-photon processes in an ideal case. In this ideal 
case, with the perfect pulse-duration and frequency-matching 
conditions, we can prepare the perfect target state. However, 
in practical cases, the system cannot avoid environmental ef¬ 
fects. Moreover, the imperfection of the parameters chosen 
also affects the fidelity of the target state. For example, dif¬ 
ferent N describe different Bessel functions for effective cou¬ 
pling strengths between the cavity field, the two-level system, 
and the classical driving field. Then the optimization for the 
target state will also be different. For concreteness, as an ex- 
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ample, let us study the effects of both the environment and 
imperfect parameters on the target state 

|^02) = ^(|0) + |2» Iff), (43) 

in the displacement picture, whose density matrix operator 
can be given as 

P 1 = |- (44) 


We also assume that the terms with the Bessel function for 
TV = — 1 are chosen for the state preparation. But other terms 
with the Bessel function order TV' ^ TV are also involved. 
Thus we have to choose uj z ^ ton/ (TV' — TV) to minimize the 
effect of these terms. Among all the terms with the Bessel 
function order TV', the dominant ones are those with TV' = 
0,1,±2. That is, the chosen u> z has to satisfy the condition 
uj z ^ nw, nu>/2, and nco/ 3. 

To study the environmental effect on the state preparation, 
we assume that the dynamical evolution of the system satisfies 
the following master equation 

p = -i[H,p]+£ g \p]+C c \p], (45) 


when the environmental effect is taken into account, where the 
Hamiltonian H is given by Eq. ([2} and 


£q [P\ 
[p] 


l r yij( 2 ViiPVji-°nP- P°jj), (46) 

l>j>2>0 

— (2 apa) — a) ap — pa)a) , (47) 


describe the dissipation of the qubit and the cavity field, re¬ 
spectively. Here p is the reduced density matrix of the qubit 
and the cavity field. And a z j = |*)(j|, where we define 
|0) = |ff), and |1) = |e). The operators < 7 ^ are given by 

Vij = R y ( 9 ) crijRy ( 9 ), 

with R y (9 ) = exp (—i0a y / 2) and 8 = arctan (uj x /aj z ). This 
is because we have used the eigenstates of a z as a basis (per¬ 
sistent current basis) to represent the Hamiltonian of the qubit. 
Note that 710 is the the relaxation rate, while 711 and 700 are 
the dephasing rates. The parameter n is the decay rate of the 
cavity field. In the following calculations, we assume 700 = 0. 

We first neglect the environmental effects and just study 
how the unwanted terms with the Bessel function for TV' ^ 
— 1 affects the fidelity for different parameters x,i = ‘2{1 ,i/oj ( i 
and 77 = 2 ff/w of the driving field and the cavity field when 
the target state in Eq. ( 1431) is prepared. We define the density 
matrix 


P D = ^ 02 ) (i’ 02 1 = D\p/2) p 1 D(j]/2) (48) 

which is the ideal target state in the original picture. The ac¬ 
tual target state generated in the original picture is denoted by 
the density operator p A when the effect of the unwanted terms 
is taken into account. The Fidelity for the target state is then 
given by 


(49) 


TABLE II: The fidelities of the target state are listed for different 
values of the parameters Xd = 2fld/uid and 77 = 2g/u>. Here we 
have chosen u z /2n = 19.5 GHz, u) x = 1.6 GHz, and u = 2 GHz. 



Lamb-Dicke parameter 77 = 

to 


0.330 

0.590 

0.850 

1.110 

1.370 


0.265 

0.240 

0.219 

0.237 

0.347 

0.332 


0.525 

0.281 

0.282 

0.621 

0.719 

0.734 

Xd — 2 

0.785 

0.445 

0.416 

0.722 

0.827 

0.841 


1.045 

0.369 

0.519 

0.780 

0.857 

0.877 


1.305 

0.335 

0.530 

0.791 

0.886 

0.879 


Let us now take the parameters oj/ 27t = 2 GHz and 
u> z /2it = 19.5 GHz as an example to show how the parame¬ 
ters affect the fidelity. The highly symmetry-broken condition 
is satisfied by taking, e.g., u) x /2tt = 1.6 GHz. We have listed 
the fidelities for different 77 and in Table HU from which we 
can find that larger values of 77 and Xd are more likely to in¬ 
duce a higher fidelity. Because in the range considered for 77 
and Xd, a larger . 17 / can enhance the desired term through mak¬ 
ing J_ 1 (xd) larger [see Fig. EJa)] while 77 achieves the same 
goal by enhancing the Rabi frequency for |0) \e) <—> |2) |ff) 
(see Fig. [7] in Appendix lAl). 

In Table HU the largest fidelity is F rn = 0.886 which occurs 
at the optimal parameters 77 = r] m = 1 . 11 , and Xd = x™ = 
1.305, where we have also obtained the total time T m = 1.82 
ns for generating the target state. From the above numer¬ 
ical calculations, we show that the fidelity of the prepared 
target state is significantly affected by the parameters of the 
qubit, cavity field, and driving field. Note that the fidelities 
in Table [TT| may not be satisfactory for practical applications 
in quantum information processing, which may require fideli¬ 
ties approaching 100%. However, the fidelity can be further 
optimized by carefully choosing suitable experimental param¬ 
eters. For instance, 77 = 1.5 and Xd = 1.305 would produce 
a more desirable fidelity of 0.9143, and it is still possible to 
obtain much higher fidelities when related parameters are fur¬ 
ther optimized. We should also mention that the effect of the 
unwanted terms can be totally avoided if for each generation 
step, the control pulses for the driving frequency aid, driving 
strength f Id, driving phase 77 , and the pulse duration t are all 
perfectly designed to compensate the effect of the unwanted 
terms. 

Now we study the environmental effect on the fidelity of 
the prepared state by taking experimentally achievable param¬ 
eters, e.g., 7 io/ 27 t = k/ 2 tt = 1 MHz and 711 = 2 MHz. We 
also choose 77 = rjm and Xd = x and other parameters (i.e., 
Wj, oj z , and ui) are kept the same as in Table HP Now the fi¬ 
delity we obtain via numerical calculations is F' n = 0.8775. 

The Wigner function represents the full information of the 
states of the cavity field and can be measured via quantum 
state tomography |64]. The Wigner function of the cavity field 
has recently been measured in circuit QED systems [65, 66 ]. 
To obtain the state of the cavity field, let us now trace out the 
qubit part of the density operator for the qubit-cavity compos- 


F = Tr{p A p D }. 
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ite system using the formula 

p p c = Tr q pP = (g\ p p \g) + (e\ p p \e), (50) 

where p refers to either D, I, or A. Here, p 1 is the ideal 
target state in the displacement picture, p D is the the ideal 
target state in the original picture, and p A is the actual tar¬ 
get state in the original picture. Therefore, p v c is the cavity 
part of the qubit-cavity-composite state p p . It should be em¬ 
phasized here that the actual state p A denotes the generated 
target state in the original picture with the same parameters 
as in the ideal case, but including the effects of both the en¬ 
vironment and unwanted terms. By definition, given an arbi¬ 
trary density operator p 1 the Wigner function VV (j3, (3*) and 
the Wigner characteristic function C w (A, A*) have the fol¬ 
lowing relations 

C w (A, A*) = Tr {pexp (Aa f — A*a)} , (51) 

W (ft / 3 *) = C w (A, A*) exp (-A/T + A*/3). 

(52) 

Moreover, if p is expanded in the Fock state space, i.e, 

p = '^2pmn\m) (n\, (53) 


then we have the Wigner function of p given by 

W(0,p*)='52p mn Wm n {P,P*), (54) 

ran 


where 


2 n-m + l m 


(-i r x %/3 n 


Wmn (j3,/3*) = 


xe -m 2 L £- m) ( 4 1 / 3 | 2 ), 

n m — n-t-1 , . / 1 „ 

2 __ f _1 \ n I n \_ o * m—n 

7T A ' V m - ^ 

_ xe- 2 ^ 2 Ll m - n) (4|^| 2 ), 


m < n, 

m > n. 

(55) 


As shown in Eq. ( l54l ). the Wigner function and the den¬ 
sity operator can in principle be derived from each other, 
which are closely related by the function W mn (f3,/3*). If 
{y^mn (/3,/3*)} are taken as the basis functions, then p mn 
can be considered as the spectrum of W (/3, /3*). Moreover, 
if we define p D = D{a) p D\a) and its Wigner function as 
W D ( j3 , f3*), through the definitions in Eq. (BTT > and Eq. (l52l) . 
we can easily obtain 


W D (f3, (3*) = W (/3 — ct, /?* — a*). (56) 


It is clear that the displacement operator D{a) displaces 
the Wigner function by a in the coordinate system. Since 
Pc = D{r] m /2) pc D^(r] m /2), the Wigner function for pf, 
W®(/3,/3*) and that for p ! c , W,[(/3,/3*) must have the rela¬ 
tion Wc (P,/3*) = W/(/3 - Vm/2, [3* - Pm/2). Therefore, 
Fig. [6ja), i.e., the figure for and Fig. |6jb), i.e., the fig¬ 
ure for VV ( /J , are in fact of the same profile except that there 


is a horizonal translation between them. In Figs. [6j a,h,c), the 
vertical dashed line that goes through the maximum value of 
the Wigner function indicates the horizonal component of its 
central position. Since the displacement operator between 
Fig. [6ja) and Fig. Ob) is I) (:}„, /2), then the amount of the 
translation is r/ m /2 = 0.555. When including the environ¬ 
ment and unwanted terms. Fig. O c ) shows how the Wigner 
function becomes different from Fig. Ob)- We can determine 
that the displacements of Fig. Ob) and Fig. 0<-') are basi¬ 
cally the same. But a careful comparison shows that the hori¬ 
zonal central position of Fig. O c ) is 0.4745, 0.0805, which 
is less than that of Fig. Ob), which is rj m /2 = 0.555. We 
think this small difference can be mainly attributed to the ef¬ 
fect of the environment and unwanted terms. FigureO c ) also 
shows local twists as well as a global rotation compared with 
Fig. Ob). The global rotation represents the average phase 
noise, while the local twists represent the corresponding fluc¬ 
tuations. Though both the environment and unwanted terms 
both affect the fidelity of the states prepared, our calcula¬ 
tions show that under the specified parameters, the role of 
the unwanted terms is dominant when the imperfect pulses 
are applied to state preparation, since the generation time 
T m = 1.82 ns is far from inducing serious decoherence at the 
specified decay rates, which is well manifested by the poor 
fidelity reduction F rn — F’ m = 0.0085. Recall that F m is the 
fidelity obtained using the optimal paramters in Table. Ullwhen 
only including the effects of the unwanted terms, while F' rn is 
the fidelity obtained using the same parameters, but with the 
effects of both the environment and unwanted terms consid¬ 
ered. 


VI. DISCUSSIONS 


Let us now discuss the feasibility of the experiments for 
the generation of nonclassical microwave states using super¬ 
conducting quits interacting with a single-mode microwave 
field. The frequency of the qubit cannot be extremely large. 
Thus the maximum photon number in multiphoton processes 
is limited by the ratio between the frequency oj z of the 

qubit and that of the cavity field w. This means that the qubits 
should be far away from the optimal point for the flux and 
charge qubits when the microwave states are generated using 
our proposed methods. This might be a problem for the prepa¬ 
ration of arbitrary superpositions, because the coherence time 
becomes short when the flux or charge qubit deviates from the 
optimal point. However, for the particular number state |n), 
there is no requirement for the coherence and thus it should 
be more efficient, because we need only to prepare the qubit 
in the excited state, and then the state \n) can be prepared 
via an n-photon red-sideband excitation. We know that the 
phase |6d, [6ll and Xmon Hi qubits are not very sensitive 
to the optimal point. Thus the proposal might be more effi¬ 
cient for these qubis coupled to a microwave cavity. It should 
be noted that the imperfect pulse can significantly affect the 
fidelity. We thus suggest that enough optimization be imple¬ 
mented to reach an acceptable fidelity. 
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FIG. 6 : Wigner functions (a) W c f (/?, /?*), (b) W<? (/3,/3*), (c) 
(13,13*) respectively for (a) pi, (b) p®, and (c) pf. Here 
x = Re(/3), and y = Im(/3). In the above figures, the param¬ 
eters are chosen as ui z /27r = 19.5 GHz, u* = 1.6 GHz, oj = 
2 GHz, 721/271 = k/2tv = 1 MHz, 711 = 0, 722 = 0.1 MHz, 
77 = g m = 1.11, and Xd = x™ = 1.305. We have used a ver¬ 
tical dashed line in (a), (b), and (c) to highlight the displacement 
of the central point, also the maximum point, of the Wigner func¬ 
tions. As shown in (a), the central point of Wi (/3,/3*) is the ori¬ 
gin. Since p c = D(g m / 2) pi (gm/2), the Wigner function 

for p ®, W,? (/3,/3*) and that for pi, V\?Z (/3,/3*) h ave the relation 
(/3, X) = VVc (/3 - r)m/2, /3* - g m / 2). This is the very rea¬ 
son why (a) and (b) exhibit the same profile as well as a horizonal 
translation. The exact value of this translation length is, of course, 
r/m/2 = 0.555. Compared to the ideal target state in the original pic¬ 
ture, i.e., Wf (0, 0*) in (b), the actual target state W(0, X) in ( c ) 
possesses nearly the same central point. But due to the effects of en¬ 
vironment and unwanted terms, there appears, in (c), a new feature 
of small local twists and a global rotation. 


VII. CONCLUSIONS 

We have proposed a method to prepare nonclassical 
microwave states via longitudinal-coupling-induced multi¬ 
photon processes when a driven symmetry-broken super¬ 
conducting qubit is coupled to a single-mode microwave 
field. With controllable /c-photon processes in a SQ with a 
symmetry-broken potential energy, only n max + 1 steps are 
needed to synthesize the superposition of Fock states with 
the largest photon number nm ax . However, in contrast to the 
method used in Refs. 15811611] . with one-photon processes in 
the SQ inside the cavity, 2n max steps are needed to synthe¬ 
size the same state. Moreover, using A;-photon processes, a 
/c-photon Fock state | k) can be generated with just two steps, 
while with one-photon processes, 2k steps are required to pro¬ 
duce the same state. Thus, the time to generate the same state 
using multiphoton processes is shorter than that using only a 
single-photon process. Therefore the fidelity should also be 
improved. In this sense, our method is more efficient than the 
one in Refs. ll58i.l6ltl . Besides, we have provided an analyti¬ 
cal solution for the total time needed to generate a target state 

lA'-Vlmax)' 

We have made a detailed analysis of the ground state when 
the system is sufficiently cooled. We find that in the highly- 
symmetry-broken and strong-coupling case, the ground state 
can still be regarded as the vacuum state in the displacement 
picture. The displacement effect on both the initial state and 
the target state has also been studied. Generally, the displace¬ 
ment will induce a variation of the photon-number distribu¬ 
tion. But in the representation of the Wigner function, its in¬ 
fluence is just a shift of the center of the Wigner function by 
the Lamb-Dicke parameter 77 = 2 g/oj between the coupling 
strength g of the cavity field to the qubit and the frequency u> 
of the cavity field. We note that the Fock state produced in the 
displacement picture is a displaced number state in the origi¬ 
nal picture. Thus, a circuit QED system with broken symme¬ 
try in the qubit potential energy can be used to easily generate 
a displaced number state. This can be used to study the bound¬ 
ary between the classical and quantum worlds II78 h81|| . 

In summary, although we find that the nonclassical pho¬ 
ton state can be more easily produced when the symmetry of 
the potential energy of the SQ is broken, this method can be 
applied to any device with longitudinal and transverse cou¬ 
plings to two-level systems. Although the Fock state can be 
produced in any symmetry-broken qubit, the superposition of 
Fock states might be easily realizable in a circuit QED system 
formed by a phase qubit and a cavity field. This is because 
phase qubits have no optimal point, and thus not sensitive to 
the working point of the external parameter. Our proposal is 
experimentally realizable with current technology. 
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Appendix A: Detailed steps for generating the nonclassical state 


if we select an N and Xd such that Jn ( Xd ) > 0. Here, ui^ d \ 
and t n are, respectively, the driving frequency, driv¬ 
ing phase, and time duration for each generation step. From 
Eq. ||A3), we know that 


If we substitute Eq. (124b . Eq. (126b . and Eq. (1321) into 
Eq. (PTT) . then the following relations can be obtained, e.g., 
for the generation of |'i/ , o)> i.e., for the step n = 0 , 


a 


(n) 

fcO 


rt()lmax) 

L 'fcO 


= \C k \, k < n. 


and hence 


C 0°J = ex P ( ia 0 e) sin ( 


U°’° 

“JV,C 


y- 1 ) 

°0 g > 

(Al) 

yr(n-l) 

°0e 

/ n—1 

= >-E|<r 

O 0 ’ 0 

s ‘JV,c 

to 

y _ t) 

°0g ) 

(A2) 


\ fe= 0 


1/2 


= E to? 


(A13) 


1/2 


v k—n 


with | Cq 9 1} | = 1 , which is determined by the initial condi¬ 
tion. However, for the generation of \ijj n ) with n > 1 from the 
state we can obtain the following relations for their 

coefficients 


y«) — pvn (in n 1 r^” 1 ) 

L 'kg ~ eX P \ la kg) °kg 
C ^ng = ex P (*<g) sin ( 
Cte = eX P (*« 0 e) C0S ( 


O n ’° 

si JV,r 


o".° 

“iV,r 


b nJ v^Oe > 

i \ y ra_ t) 

b n ) ^0e i 


(A3) 

(A4) 

(A5) 


with k < n — 1. Here, the phases and for n = 0 are 
determined by 


a (0) = 

a 0g - 


•Ed . / , 

—sm [u K d 't 0 + 


,(o). 


^ 0) ) - y sin (<^ 0) ) 


„( 0 ) _ x d. (' ( 0 ), ,( 0 )\ Xd 

a 0e ~ ^~ sm Y d to + <Pd ) -TF Sm 

, 7T 
+ <^ 0 - 2 - 


+ 


) 


The other phases with forn > 1 are given by 

v (n) _ f (n) 


a kg = 


ysin (+ 4 " } ) - ysin (^" } ) 


- kuit n , 

ang 1 = ysin + 4'° 

71 

TlUjt n (Pn “ 


(n) 

a 0e = - 


Xd . ( 

T sm l 1 


2 ’ 

, ,(")+ _i_ ,/>(") 

OJ d t n + cp d 


Xd . 
—sm 
2 


Xd . 
-sm 


) 


Uztp 

2 ’ 
(A 6 ) 

_ Uztp 
2 

(A7) 


0J z tn 

2 

(A 8 ) 

OJztn 

2 

(A9) 

u z t n 

2 ’ 

(A10) 


where A; < n — 1. In Eqs. (1A1HA5I) . 

arg (^w°c) = - *3 n = 0 , 

(^ 7 V°r) = - /V< / > i n) - {n + 1) 7T, n > 1, 


arg | 


if we select an At and x d such that Jjv (x d ) < 0 , and 


( o 0,0 \ a r 

arg ( c ) =N$y, 


n = 0 , 


arg ( ) - rwr, n > 1, 


(All) 


(A12) 


(A14) 

Then, from Eqs. (1A 1 1lA2b and Eqs. (IA41 lA5b . we respectively 
have 


n°’° 

“JV,c 


o n -° 

si JV,r 


to = arccos 


= arccos 


t.„ = arcsin 


= arcsm 



(A15) 


+ 2ln. (A16) 


where l is an arbitrary integer. Using Eqs. (ED E3 and 
Eqs. (1A4IIA5I) . we derive 

arg (^ 0 ° } ) - arg (c^ = a ( 0 °J - a ( 0 °J + 2ln, (A17) 
arg (C^) - arg ( C^ = a$ - a^ + 2ln. (A18) 

If n = rimax in Eq. (IA18b . then C^"' 1 = 0, with no definition 
of the phase angle. Thus we can assume that arg = 0, 

(n) (n) 

without affecting the final result. Here, Cng and C (lf , can be 
obtained through the following recursion relations 


y™-t) _ 

°0e ~ 


C 


(n) 


exp(ia^) cos(|n"’“ |t„) 


, 1 < n < n n 


c ( ”) 


exp ) sin( | | t n ) 




r<t n ) 

kg 


exp 


(»£’) ’ 


k < n — 1 , 


y _ i) _ 
°o g - 


C 


(o) 

Oe 


exp sin ( U^° c f 0 ) 


(A19) 

(A20) 

(A21) 


In Eq. ( 1A 19b . distinguishing the case when n = n lnax from the 
other ones is needed to avoid the apperance of 0/0. Though 
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FIG. 7: (color online) Reduced Rabi frequency from 

Eq. <A22l . as a function of 77 , for n = 1,2,3,4, and 5, respectively. 
The square on each plot denotes the point that achieves the largest 
Recall that V = 2g/uj, where g is the qubit-cavity coupling 
constant and oj is the frequency of the single-mode cavity field. Thus, 
77 is the normalized coupling. 



FIG. 8 : (color online) Plot of log yTj as a function of 
the Lamb-Dicke parameter 77 = 2g/uj for |f(Vi max ) = 
5 ^n=o x l n ) ® I®) /\/ n max + T, with different maximum photon 
number n max . Recall that T = T\u>Jn\ /2 is a normalized time 
and T is the total time to generate the desired target state. The star 
on each curve shows the optimal point where the normalized gener¬ 
ation time reaches its minimum. 


Eq. (1A2 1 b implies that Cg g may have a definite phase, such 
a phase could only add a global phase factor to the target state. 
So it is convenient to directly specify Cq ~ 11 = 1. 

We define the reduced Rabi frequency as 


O 


n, 0 
N,/3 


2 


n 


n, 0 
N,P 


I^xJn I 



(A22) 


in order to study its dependence on 77 . From Eq. ( IA22L we 
can obtain the optimal Lamb-Dicke parameter 

t?n,o — (A23) 


that achieves the largest reduced Rabi frequency 




n, 0 
N,/3, o 




(A24) 


which is also the point that makes 
illustrated in Fig. [7] We can also verify 


qO,u 


q0,71— 1 

lL N,p 


as 


lim 

n—t 00 


qti+1,0 


o".° 

N,0,o 


= lim 



n + 1 


= L 


(A25) 


with limn^oo {rjn+i,o/iln, 0 } = 1- This means that when 
the photon number n increases, the optimal points for the 
Rabi frequencies between the zero-photon state and different 
n photon states tend to approach each other infinitesimally. 
But for low photon numbers the optimal points are still distin¬ 
guishable from each other. 

Let us calculate the total time T for generating the target 
state 


Umax 

T = ^2 t n 


n =0 


:- 7 —r arccos (| C 0 \) exp 

l^x^N I 



Umax 

+ £ 


n =1 


2V^- 

\u x Jn \ r] n 


arcsin 




(A26) 


where 77 = 2g/u is the Lamb-Dicke parameter and we have 
omitted excessive cycle periods for each step. By taking the 
derivative of T with respect to 77 , we can find all the extreme 
points of 77 , which satisfy the following equation 

^max4' 1 

A„77 n ““ +1 -" = 0. (A27) 

n=—1 

The coefficient A n has been given in Eqs. ( lA31l )-( lA34l ). 
Further selection among these extreme points and the 
experimentally-constrained boundaries of 77 can yield the op¬ 
timal Lamb-Dicke parameter 77 opt , which will lead to the least 
generation time T opt . Once T opt is reached, in principle, the 
influence of the environment on the target state fidelity will be 
minimized. 























































14 


Similarly to 


Q°’ ra 


, we define 

T{ri) =T\u x J n \/2. 


(A28) 


The curves of log (t ( 77 )^ which for particular states have 

been plotted in Fig. [8] with a star on each curve to label the 
point where the generation time reaches its least value. 

The normalized time needed to generate a target state is 


T = ^ t n = arccos (|Co|) exp (rj 2 /2) 

| C n \ \ e ^ 2 / 2 


n —0 


+ ^2 


arcsm 


,\jTr k ~n\Ck\ 2 ) vr 


whose extreme points still satisfy 

Umax “ 1 “ 1 


E A » 


V 


i^max'1‘1 ^ — 


= 0 , 


n=— 1 

for unbound rj, where if n max = 0, 


— 


if rtmax — I ■ 


A n — 


arccos (|Co|), n = — 1, 
0, others, 


arccos (|Co|), n = — 1, 

P n + 1 , n = 0, 

- (n - 1) Pn- 1 , n = 2, 

0, others, 


(A29) 


(A30) 


(A31) 


(A32) 


if n max — 2, 


A-n — < 


arccos (|Co|), n = — 1, 

Pn+i, 0 < n < 1, 

- (n - 1 ) -P n -i, 2 < n < 3, 

0, others, 


(A33) 


and for other cases, we have 


arccos (|Co|), n = — 1, 

P„+i, 0 < n < 1, 

Ar). < Cn+i (rt, 1) P n — 1 , 2 E rr E ^max 1? 

(rt 1) Pn—lj ri max E ri E ri max -{- 1, 

0, others. 


(A34) 


Here we have used the abbreviation 


ICE 


P„ = V it! arcsin 


/ V^ n max I/^T I 2 

, V Z-^k—n \^k\ 


(A35) 
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